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Abstract. In artificial neural networks, neurons are usually implemented with 
highly dissipative CMOS-based operational amplifiers. A more energy-efficient 
implementation is a “spin-neuron” realized with a magneto-tunneling junction (MTJ) 
that is switched with a spin-polarized current (representing weighted sum of input 
currents) that either delivers a spin transfer torque or induces domain wall motion in 
the soft layer of the MTJ. Here, we propose and analyze a different type of spin-neuron 
in which the soft layer of the MTJ is switched with mechanical strain generated by a 
voltage (representing weighted sum of input voltages) and term it straintronic spin- 
neuron. It dissipates orders of magnitude less energy in threshold operations than the 
traditional current-driven spin neuron at 0 K temperature and may even be faster. We 
have also studied the room-temperature firing behaviors of both types of spin neurons 
and find that thermal noise degrades the performance of both types, but the current- 
driven type is degraded much more than the straintronic type if both are optimized 
for maximum energy-efficiency. On the other hand, if both are designed to have the 
same level of thermal degradation, then the current-driven version will dissipate orders 
of magnitude more energy than the straintronic version. Thus, the straintronic spin- 
neuron is superior to current-driven spin neurons. 


Keywords : spin-neuron, neural network, straintronics, nanomagnets 

1. Introduction 

The building blocks of neural computing architectures are ‘neurons’ usually connected 
to each other and to external stimuli through programmable ‘synapses’. The transfer 
function of a neuron can be expressed as 



( 1 ) 


where / is some nonlinear function, uy-s are programmable weights of synapses, Xj-s 
are the input signals (representing dendrites), b is a fixed bias and O is the output 
(representing a neuron’s axon). In threshold operations, / mimics a Heaviside unit step 
function whose value is 1 if the argument is positive and 0 otherwise. 
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Neural networks have myriad topologies, such as cellular neural network [1], feed¬ 
forward network [2], convolutional neural network [3], hierarchical temporal memory 
[4], etc. However, the basic unit of computing i.e., the neuron, remains more or 
less invariant across all topologies and its operation is governed by Equation (1). In 
a conventional neural network, CMOS operational amplifiers carry out the threshold 
operation of Equation (1) [5] and dissipate exorbitant amounts of energy. To a large 
extent, this has stymied the progress of neural computing. Alternate implementations 
to lower the energy dissipation have been proposed in recent years [6, 7, 8] and utilize a 
magneto-tunneling junction (MTJ) whose soft layer is an anisotropic nanomagnet with 
two stable magnetization orientations. Input variables are encoded in spin-polarized 
currents that are summed with variable weights to produce a net spin polarized current 
which is driven through the nanomagnet. When the net current exceeds a threshold 
value, the magnetization of the soft layer rotates from one stable orientation to the other, 
thereby changing the resistance of the MTJ abruptly. This implements the threshold 
bring behavior of a neuron. These types of artificial neurons have been termed ‘spin- 
neurons’ and unlike CMOS-based neurons, they are ‘non-volatile’ since the final state of 
the neuron can be stored in the magnetization state of the nanomagnet (and therefore 
the resistance of the MTJ) after the device is powered down. In Ref. [8], the soft layer 
of the MTJ is a nanomagnet possessing perpendicular magnetic anisotropy (see Figure 
3(a) of Ref. [8]) and it is switched with spin-polarized current generated via the giant 
spin Hall effect [9] which induces domain wall motion. This type of spin neurons belongs 
to the general class of (spin-polarized) current driven artificial neurons. 

In this paper we propose and analyze a different type of spin-neuron implemented 
with MTJs having soft layers that are magnetostrictive or multiferroic nanomagnets 
and whose magnetizations are flipped with mechanical stress/strain generated by a 
voltage. We call them ‘straintronic spin-neurons’ and they are voltage-driven as opposed 
to current-driven. This has the advantage of further reducing the energy dissipation 
during bring. Switching of multiferroic nanomagnets with voltage-generated stress 
has been proposed and/or demonstrated by many groups [10, 11, 12, 13, 14] and is 
particularly useful for writing bits in non-volatile memory [15, 16, 17, 18, 19, 20]. It 
can be also harnessed for logic applications [12, 21, 22, 23] and results in exceptionally 
low dissipation. Here, we propose it for neural applications. We compare the energy- 
efficiency of a straintronic spin neuron with that of a traditional current-driven spin 
neuron and show that the former is more energy efficient. Finally, since magnetization 
dynamics is vulnerable to thermal noise, we study the operation of spin-neurons at room 
temperature in the presence of thermal noise and compare that with 0 K operation to 
assess the degree of thermal degradation. As expected, thermal noise has a deleterious 
ebect on the threshold behavior and seriously degrades the abruptness of the bring 
action. The degradation is far worse for the current-driven type than for the straintronic 
type. 
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Figure 1 shows the schematic of a straintronic spin-neuron with programmable synapses. 
Inputs Xj-s and the bias b are in the form of voltages V)-s and b, the latter being 
realized with a constant current source I [b = I (Ri || R 2 || r\ || r 2 || • • • || r^-i || tn)]- 
The voltage appearing at node P is dropped across the piezoelectric layer underneath 
the (shorted) contact pads A and A'. It is a weighted sum of input voltages and bias, 
and is given by 

N 

Vp = 5>iVi + &, (2) 

2—1 

where w, = „J R ,f ^^-T } r ’A, 1 ^7 ^ . The resistances R\,R 2 are the resistances of 
the piezoelectric layer underneath the contact pads and r,-s are the series resistances 
(connected to the input terminals) that implement the programmable weights. Equation 
(2) is obtained from voltage superposition. 

The magneto-tunneling junction (MTJ) in Figure 1 is the central unit of the neuron. 
It has a hard nanomagnet, a spacer layer and a soft magnetostrictive nanomagnet in 
contact with the piezoelectric. All nanomagnets are shaped like elliptical disks. A 
bias magnetic held B in the plane of the soft nanomagnet directed along its minor 
axis makes the magnetization orientation of the soft nanomagnet bistable, with the two 
stable directions shown as Ty and Tx which subtend an angle of ~90° between them 
[15]. The hard nanomagnet is implemented with a synthetic anti-ferromagnet and its 
two stable magnetization orientations are roughly along its major axis because of the 
extremely high shape anisotropy that this nanomagnet possesses. The hard nanomagnet 
is placed such that its major axis is colliuear with one of the stable magnetization 
orientations of the soft nanomagnet (say Ty), resulting in a “skewed MTJ stack” where 
the major axes of the two nanomagnets are at an angle. The hard nanomagnet is then 
magnetized permanently in the direction that is anti-parallel to T y. Thus, when the soft 
nanomagnet is in the stable state Ty, the magnetizations of the hard and soft layers of 
the MTJ are mutually anti-parallel, resulting in high MTJ resistance, while when the 
soft nanomagnet is in the other stable state Tx, the magnetizations of the two layers are 
roughly perpendicular to each other, resulting in lower MTJ resistance. The ratio of the 
high-to-low MTJ resistances is approximately 1/ (1 — 771 ^ 2 ), where the r/-s are the spin 
injection/detection efficiencies at the interfaces of the spacer with the two nanomagnets. 
We assume that at room temperature, ip = rj 2 ~ 70% [24] and therefore the resistance 
ratio will be roughly 2:1. Higher resistance ratios exceeding 6:1 have been demonstrated 
at room temperature [25], but we will be conservative and assume the ratio to be 2:1. 

The electrodes A and A' are placed on the piezoelectric layer such that the line 
joining their centers is parallel to Ty and hence also to the major axis of the hard 
nanomagnet. Their lateral dimensions are of the same order as and the inter-electrode 
separation is 1-2 times the PZT thin him thickness. When voltages are applied between 
the electrode pair and the (grounded) conducting substrate, electric fields are generated 
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Using superposition 

N 

V p = '£ J w i V t + b 


where w, = 


■R; ll-S 2 Ik, il r, ||•»«il11^, ||••«II 


mdb = I(R, !!«, ||r,|| r 2 ll***!^ II***!!»v) 


K = i s R m = fh'p) = f\ 5 >>+ b I = /[ X>T + b 

=>0 = /f 


Figure 1. Schematic of a straintronic spin-neuron implementing a step transfer 
function. The artificial synapses are realized with the passive resistors r± ■ ■ ■ r n . 


underneath the electrode pads in the PZT layer as shown in Figure 1. They produce 
out-of-plane compressive strain and in-plane tensile strain or vice versa, depending on 
the polarity of the voltage in the PZT layer below the electrodes [26]. These strain fields 
interact and produce biaxial strain between the electrodes (tensile along the line joining 
the electrodes and compressive in the perpendicular direction, or vice versa) [26], which 
is almost completely transferred to the magnetostrictive soft magnet since the latter’s 
thickness is much smaller than that of the strained PZT layer. 

If the magnetostriction coefficient of the soft magnet material is positive (e.g. 
Terfenol-D), then sufficient compressive stress resulting along the line joining the 
electrode centers, i.e. in the direction of Ty, will rotate the soft layer’s magnetization 
to Tp, while sufficient tensile stress will keep the magnetization aligned along 4/y. The 
situation will be opposite if the magnetostriction coefficient of the soft magnet is negative 
(e.g. cobalt). Since, for either sign of the magnetostriction coefficient, the sign of the 
stress (compressive or tensile) depends on the polarity of the voltage applied between 
the electrodes and the grounded substrate, the magnetization of the soft magnet can be 
aligned along either of the two stable orientations Ty and Tp at will by simply choosing 
the voltage polarities (of the inputs and bias voltages). 

There is a minimum stress (compressive/tensile) required to switch the 
magnetization of the soft nanomagnet of the MTJ from one orientation (say, Ty) to 
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the other (say, 'hjJ because the two stable states are separated by an energy barrier 
[15, 18, 19, 23] that needs to be overcome by stress to make the switching occur. At 
0 K temperature, this feature gives rise to a sharp threshold in the switching behavior 
and makes it possible to mimic the sudden bring behavior of a neuron. The energy 
barrier (and hence the threshold stress) depends on the permanent magnetic held B 
and the size and shape of the soft nano magnet, if we ignore dipole coupling with the 
hard nanomagnet and any magneto-crystalline anisotropy (the magnets are assumed 
to be amorphous). These parameters determine the effective in-plane energy barrier 
between the two stable magnetization states T y and Jw that must be overcome by 
stress to switch the magnetization of the soft nanomagnet from one state to the other 
and therefore the MTJ resistance. The minimum stress needed for switching (also called 
the ‘critical stress’) at 0 K can be found by equating the stress anisotropy energy to the 
effective in-plane energy barrier. 

The critical stress gives rise to a critical voltage V c for switching at 0 K. When 
the total voltage Vp, appearing at node P, due to all weighted inputs and the bias, 
exceeds the critical voltage, the MTJ resistance switches abruptly because the soft layer’s 
magnetization rotates from one stable orientation to the other. If we bias the MTJ with 
a constant current source Ip as shown in Fig. 1, then the voltage across the MTJ is 
Vmtj = IbRmtj where Rmtj is the MTJ resistance that has a non-linear dependence 
on Vp (abrupt switching when Vp exceeds V c ). Because there is virtually no electric field 
in the PZT layer directly underneath the magnet, we can ignore any potential drop in 
this region and write the output voltage in Fig. 1 as Vo ~ Vmtj = IbRmtj which has a 
non-linear dependence on Vp and hence exhibits a threshold behavior. Using Equation 
(2), we can now write 

WiV + b'j , (3) 

which replicates the neural behavior. 

3. Simulation of neuron firing with the Landau-Lifshitz-Gilbert (LLG) 
equation 

We will design our straintronic spin-neuron by choosing the dimensions of the 
nanomagnet and the magnitude of the permanent magnetic field to select the critical 
voltage V c , and therefore the firing threshold. We choose a soft nanomagnet of major 
axis a = 100 nm, minor axis b = 42 nm and thickness d = 16.5 nm, which ensures that it 
has a single ferromagnetic domain [27]. When stress is applied on the soft nanomagnet 
by applying the voltage Vp at the electrodes A and A', which is the weighted sum of the 
inputs and bias, the magnetization vector of the soft nanomagnet experiences a torque 
that makes it rotate and ultimately switch the MTJ resistance. The torque depends on 
the shape anisotropy energy of the soft nanomagnet, the permanent magnetic field and 
the stress. 


N 




=/(u)=/nr 


, i =1 
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The shape anisotropy energy is a function of the time-dependent magnetization 
orientation during switching and can be written as 

E sh (t) = E s i(t)sm 2 6' (t) + E s2 (t)sm20’(t) + M 2 (N d _ yy + N d _ zz ) 

E si(t) = (y) VM 2 (2 N d _ xx - N d _ yy - N d _ zz ) cosV(f) 

E s2 (t) = (^) QM 2 (N d _ zz - N d _ yy ) sin 0'(t), (4) 

where 6'{t) and <p'{t) are, respectively, the instantaneous polar and azimuthal angles 
of the soft nanomagnet’s magnetization vector in the primed reference frame shown in 
Figure 1. The unprimed reference frame is such that the z-axis coincides with the soft 
nano magnet’s easy axis and y-axis with the hard axis. The primed reference frame 
is obtained by rotating anticlockwise about the x-axis by 45°. The quantity M s is 
the saturation magnetization of the soft nanomagnet, N d _ xx , N d _ yy and N d _ zz are the 
demagnetization factors that can be evaluated from the nanomagnet’s dimensions [28], 
Po is the permeability of free space, and hi = (7t/4 )abd is the soft nanomagnet’s volume. 

The permanent magnetic field contributes an additional term to the potential 
energy of the soft nanomagnet given by 

E m (t) = —=M s VlB [cosO'(t) — sin 0'(t)sm (j)'(t)]. (5) 

v 2 

The locations of the total potential energy ( E sh (6(t ), </>(£)) + E m (0{t) 1 </>(£))) minima 
(in the absence of stress) determine the stable magnetization orientations of the soft 
nanomagnet and are found by minimizing the total potential energy with respect to 
0(t),(p(t). In our case, the two stable orientations (the two degenerate energy minima) 
turn out to be T/ 1 | (6 — 0\ — 46.9°, 0 = <f>i — 90°) and Tx (0 = 0 o = 133°, 0 = (fro = 
90° = 0i) with an angular separation 7 of 86.1° (Figure 1) between them when B is 
0.14 T. The in-plane energy barrier separating these two energy minima is 73.1 kT at 
room temperature resulting in a tiny probability (e w3,1 ~ 10 ~ 32 ) of the neuron firing 
spontaneously at room temperature by switching between the stable states 4/|| and Tx 
per attempt. With an attempt frequency of 10 15 Hz, the mean time between successive 
spontaneous firing (switching of the MTJ) would then be 10 ~ 15 x 10 32 = 10 1 ' seconds 
= 3x10' centuries. 

When voltages are applied at the electrodes A and A', the resulting stress produced 
in the soft nanomagnet contributes a stress anisotropy energy to the soft nanomagnet’s 
potential energy. Although the generated strain is biaxial, we will approximate it as 
uniaxial strain to somewhat compensate for the fact that not 100 % of the strain in 
the PZT will be transferred to the soft nanomagnet. With this assumption, the stress 
anisotropy energy is written as 

E str (t) = -^A s e(t)FHcos 2 6 |/ (t), (6) 

where A s is the magnetostriction coefficient of the soft nanomagnet, Y is the Young’s 
modulus, and e(t) is the strain generated by the applied voltage Vp(t) at the instant of 
time t. 
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The total potential energy of a stressed nanomagnet at any instant of time t is 
therefore 


E{t) = E (0'(t), <j>\t)) = E sh (t ) + E m (t) + E str (t). 


( 7 ) 


We follow the standard procedure to derive the time evolution of the polar and 
azimuthal angles of the magnetization vector of the soft nanomagnet in the rotated 
coordinate frame under the actions of the torques due to shape anisotropy, stress 
anisotropy and magnetic field by solving the Landau-Lifshitz-Gilbert (LLG) equation: 


dm(t) 

dt 


— a 


m(i) x 


dm(t) 

dt 


-hi 

ia 0 M s Q 


Ns ( 7 ) j 


( 8 ) 


where a is the Gilbert damping coefficient which depends on the soft nano magnet’s 
material, 7 is the gyromagnetic ratio (a universal constant) and t ss is the total torque 
acting on the magnetization vector and is given by 


Tss(t) 


, , / dE(t) ~ 1 dE(t)~\ 

x \dO'(t) G + smO'(t) d^it)^) 
{E^t) sin O'{t) + E^ 2 it) cos 6\t) 

—=M S VLB cos <j)'{t)}0 
v 2 

{E s i(t) sin 20'(t) + 2E s2 (t) cos 2 6\t) 


—-=M s Q,B(sm (j/(t) cos 0'{t) + sin0'(t)) 
v 2 

(3/2)A s e{t)YVt sin 2 0'(t)}<p, 


(9) 


where m(t) is the normalized magnetization vector, quantities with carets are unit 
vectors in the original frame of reference, and 

E^iit) = ^-M^Q{(N d _ yy + N d _ zz )sm2(j)'(t) - 2N d _ xx sin2(j)'(t)} 

E^it) = ^M 2 s n (. N d _ zz - N d _ yy ) cos 

At non-zero temperatures, there is an additional torque acting on the magnetization 
vector owing to thermal noise. The procedure for finding this torque has been described 
in Ref. [29] and is not repeated here. 

Solution of the LLG equation [Equation ( 8 )] yields the magnetization of the 
stressed soft nanomagnet as a function of time and steady state is achieved when the 
magnetization no longer changes appreciably with time. This yields the switching time 
(time elapsed before reaching steady-state) and the energy dissipation for any given 
stress. The procedure for finding these quantities in the presence of thermal fluctuations 
requires solving the stochastic LLG equation and is described in Ref. [29] and [30]. 

We assume that the soft nanomagnet is made of Terfenol-D which has the following 
material parameters: saturation magnetization M s = 8 x 10 5 A/m, magnetostriction 
coefficient X s = 60 x 10 ~ 5 , Young’s modulus Y = 80 GPa and Gilbert damping 
coefficient a = 0.1 [31, 32, 33]. 









Nanotechnology (2015) 


We solve the Landau-Lifshitz-Gilbert equation [Equation (8)] at 0 K (and its 
stochastic version at 300 K) to find the steady-state orientation of the magnetization 
vector of the soft nanomagnet as a function of stress in the nanomagnet and therefore 
as a function of the sum-total voltage applied at the electrode pairs. The stress is varied 
between -50 and +50 MPa. By following the procedure in Ref. [26], we compute that an 
electric held of 37.5 kV/m is required to produced a stress of 1 MPa in the PZT, which 
we assume is fully transferred to the soft magnet. Therefore, assuming that stress is 
linearly proportional to the voltage, the corresponding voltages for ±50 MPa are ±187.5 
mV if the PZT film’s thickness is 100 nm. 

At the critical stress value (or critical voltage), the soft layer’s magnetization vector 
changes abruptly from its initial stable orientation to the other causing the angle between 
the magnetizations of the soft and hard layer to change from /3l = 180° — 86.1° = 93.9° 
to Ph = 180°, or vice versa. This will cause the S-MTJ resistance to change by a factor 
of (1 ± r}iT}2Cos(93.9°)) / (1 — 771772 ) ~ 1.9 if we assume r]i = r] 2 = 0.7. 

Figure 2 shows the ratio R(V)/Rl as a function of voltage V applied at the contact 
pads, where R(V) is the MTJ resistance at a voltage V and Rl is the MTJ resistance in 
the low-resistance state. If the MTJ is initially in the high resistance state, a compressive 
stress (positive voltage) is required to drive it to the low resistance state, whereas if it 
is initially in the low resistance state, a tensile stress (negative voltage) is required to 
drive it to the high resistance state because we have assumed the soft magnet material 
to be Terfenol-D which has positive magnetostriction. Thus, the critical voltage to 
switch from high-to-low resistance is ±93.15 mV and the critical voltage to switch from 
low-to-high resistance is -93.15 mV, which produces the appearance of a ‘hysteresis’ 
in the characteristic in Fig. 2. This hysteretic behavior indicates that the device can 
also be used as a non-volatile memory. Note that the transitions are not completely 
abrupt even at 0 K temperature because sub-critical stress that is slightly lower than 
the critical stress can cause some rotation of the soft magnet’s magnetization vector 
and hence change the MTJ resistance perceptibly. This is the reason for the ‘rounded 
corners’. 

The energy dissipated during a firing event is the sum of the CV 2 dissipation 
associated with charging the capacitance of the electrodes A and A', and the internal 
dissipation in the soft magnet due to Gilbert damping. We neglect any dissipation due 
to the current sources assuming that the currents are small and the current sources are 
turned on only when the inputs arrive. Here, C is the capacitance of the electrodes 
and V is the voltage at the electrodes. The capacitance is determined from the areas of 
the electrodes, the dielectric constant of PZT and the PZT him thickness. It is found 
to be 0.88fF. Therefore, the capacitance charging dissipation will be 2 x (1/2 )CV 2 = 
0.09315 2 x 0.88fF = 7.63 aJ (there are two electrodes A and A'). The dissipation due to 
Gilbert damping at 0 K was found to be 1.2 aJ from the solution of the LLG equation as 
described in Ref. [29] and [34]. Therefore, the total dissipation during a bring event is 
8.83 aJ. This is ~450 times lower than that reported (4 fj) in Ref. [8] for current-driven 
spin neurons when the input voltage level was 100 mV and ~45 times lower than that 
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Voltage (mV) 

Figure 2. The transfer function (or firing behavior) for a straintronic spin neuron at 
0 K. When the voltage appearing at node P (in Figure 1) due to weighted inputs and 
bias is V, the resistance of the MTJ is R(V). The low resistance is Rl- 


(0.4 fj) when the input voltage level was 10 mV (our input voltage level is ~ 93 mV). 
Note that the synapse resistances rq • • • can be set arbitrarily high and hence the 
dissipation in them can be neglected. The switching (bring) delay is found to be 1.37 
ns at the threshold voltage level. 

4. Current-driven Spin Neuron Based on Spin Transfer Torque (STT) 
switched MTJ 

In this section, we discuss a spin-neuron implemented with the same type of MTJ as 
above, except this time the soft magnet is not switched with strain, but with a spin 
polarized current delivering a spin transfer torque (STT). The inputs are therefore not 
voltages, but currents. We choose the material to be CoFeB with in-plane anisotropy 
which has a Gilbert damping constant of ~0.004 [35], much lower than that of Terfenol- 
D. This material has in-plane anisotropy if the magnet thickness exceeds a few nm [36]. 
STT-driven spin neurons built with magnets exhibiting perpendicular anisotropy have 
been studied in Refs. [8] and [37]. 

Note that we chose two different materials - Terfenol-D for the straintronic neuron 
and CoFeB for the current driven neuron - because we wish to optimize both for 
minimal energy dissipation and then make a fair comparison. Terfenol-D has a very 
large magnetostriction coefficient and is hence beneficial for a straintronic neuron. A 
current driven neuron does not benefit from large magnetostriction, but benefits from 
small Gilbert damping, which is why we chose CoFeB for it. 

Figure 3 shows the schematic of a STT-based current-driven spin neuron consisting 
of an MTJ stack that is not skewed. Both magnets are elliptical in shape and each has 
two stable magnetization orientations along the major axis of the ellipse. Their major 
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Figure 3. Schematic of a current-driven spin-neuron 

axes are collinear. 

The hard magnet is permanently magnetized along one of its stable orientations. 
The soft magnet’s magnetization can be either parallel (low-resistance state) or anti¬ 
parallel (high-resistance state) to that of the hard magnet. The high-to-low resistance 
ratio in this case is (1 + r]ir] 2 ) / (1 — 771772 ) = 2.9, assuming once again that r/i = rj 2 = 0.7 
[24], This corresponds to a tunneling magnetoresistance ratio, or TMR, of ~200%. 

A negative potential applied between the hard and soft layers will align the 
magnetization of the soft layer parallel to that of the hard layer (low MTJ resistance), 
while a positive potential will make it anti-parallel (high MTJ resistance). These 
potentials cause a spin polarized current to be injected into or extracted from the soft 
layer which brings about the magnetization switching by exerting a spin-transfer-torque 
on the magnetization vector [38, 39, 40, 41]. For consistency, the soft layer’s dimensions 
are chosen such that the in-plane shape anisotropy energy barrier is 71.7 kT, close to 
that of the soft magnet in the straintronic spin neuron. This is ensured by choosing the 
major axis a = 50 nm, minor axis b = 32 nm and thickness d = 5 nm [28]. 

In the absence of spin polarized current, the soft nano magnet has only shape 
anisotropy energy and its potential energy at any instant of time t is given by, 

E sh ape{t ) = E ss (t)sm 2 6(t) + Y^ M sN d -zz 

E„{t) = (y) QM 2 {N d _ xx cos 2 cj)(t) + Nd- yy sm 2 <j){t) - N d _ zz }, (10) 

where 6{t) and 0(f) are again, respectively, the instantaneous polar and azimuthal angles 
of the magnetization vector, and M s is the saturation magnetization which is ~10.4 x 
10 5 A/m [42] for amorphous CoFeB. 

The torque on the magnetization vector at any time t due to shape anisotropy can 
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= Efait) sin 0(t)6 — E ss (t ) sin 2 6(t)<p, 


( 11 ) 


where m(i) is the normalized magnetization vector and 
E (Ps (t) = {fM] ( N d _ yy - N d _ xx )}Q sin20(f). 


Passage of the spin-polarized current I s through the nanomagnet generates a spin 
transfer torque (STT) on the magnetization vector given by [43] 


T„t(i) = s[bsin(( - 0(t))rj> - csin« - 6(t))8], 


( 12 ) 


where s = (h/2e)xh is the spin angular deposition per unit time, x is the spin 
polarization of the current, b and c are coefficients of the out-of-plane and in-plane 
components of the spin-transfer-torque. We assume x = 70% (assuming 70% spin 
injection efficiency), and b and c are 0.3 and 1, respectively. The current is passed 
perpendicular to the plane of the magnet as shown in Figure 3. The quantity ( is the 
angle subtended by the direction of spin polarization with the z-axis and it is either 0° 


or 180°. 


The Landau-Lifshitz-Gilbert (LLG) equation and its stochastic version are solved 
again to extract the STT-induced magnetization switching behavior of the soft 
nanomagnet in the absence (0 K) and presence (300 K) of thermal noise. 



(13) 
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Figure 4. The transfer characteristic (or firing behavior) for a STT-based current- 
driven spin neuron at 0 K. R(I ) is the resistance of the MTJ when the total current 
injected into the soft layer (due to weighted inputs and bias) is I and Rl is again the 
low resistance of the MTJ. 


First, we determine the current required to switch the soft layer (from the LLG 
simulations) and find it to be 127.13 /j,A (the corresponding current density is 10.11 
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MA/cnr) for a current pulse of duration 10 ns. Ref. [44] also assumed a current pulse 
duration of 10 ns and found the switching current density to be 4 MA/cm 2 for a barrier 
height of 42 kT in perpendicular magnetic anisotropy (our barrier height is 71.7 kT). 
Ref. [45] considered a CoFeB soft nanomagnet with in-plane anisotropy barrier of ~57 
kT (calculated from the reported dimensions of 170 nm x 60 nm x 1.4 nm using Ref. 
[28]) and found the switching current density to be ~ 5.5 MA/cm 2 for a switching time 
of 10 ns. 

Figure 4 shows the 0 K temperature transfer characteristic (or bring behavior) of 
the MTJ when the total input current / has been varied between -200 /iA and +200 
/iA. Note that rounded corners are more “square” here because a subthreshold current 
(unlike subthreshold voltage) hardly rotates the magnetization of the soft magnet and 
hence does not alter the MTJ resistance perceptibly. 

If we consider a typical resistance-area product (RA) of 2.1 Q-/im 2 for the MTJ 
in the high resistance state [46] and a TMR ratio of 200%, then the high-state MTJ 
resistance for our chosen dimensions becomes 1671 ohms and the low-state MTJ 
resistance 557 ohms. Energy dissipation Ep due to the current passing through MTJ’s 
tri-layered structure is given by [47] 


E r = 


dtr 


Rp + ( Rap — Rp) 


cos 0(t) 


(14) 


where Rp and Rap are the MTJ resistances in the parallel (low) and anti-parallel (high) 
state, and 6{t) is the angle between the magnetization states of the soft layer and the 
hard layer at time t. This dissipation turns out to be ~0.26 pJ. The LLG simulation 
showed that the neuron takes 14 ns to switch. The dissipation due to Gilbert damping 
in the soft magnet is a mere ~0.48 aJ, which is negligible. 

The total dissipation reported in Ref. [8] for current-driven spin neurons that 
use magnets with perpendicular anisotropy is 0.4 fj. We used magnets with in-plane 
anisotropy. There are at least two reasons why Ref. [8] could have reported a 650 times 
lower dissipation compared to what we found (0.26 pJ). First, their critical current 
density was 4 MA/cm 2 which is 2.5 times less than ours. Presumably, the lower current 
density is due to the fact that the energy barrier between stable magnetization states 
in their device might have been only ~40 kT which is what we estimate following the 
procedure in Refs. [48] and [36]. Ours was 71.7 kT. The critical current density scales 
with the energy barrier height; for example, Ref. [49] reported a critical current density 
of 8.7 MA/cm 2 for a barrier height of 67 kT. Additionally, in Ref. [8], the soft layer 
thickness was only 2 nm to maintain perpendicular anisotropy, whereas our thickness was 
5 nm and the critical current density increases with the soft layer thickness [42], These 
two factors increased the current density in our case and caused a higher dissipation. 
Second, and more importantly, Ref. [8] utilized the spin Hall effect to inject/extract 
spin polarized current from the magnets which would allow passing the charge current 
parallel to the heterointerface between the spacer and the magnets. This would allow the 
current to avoid going through the highly resistive spacer and decrease the resistance- 
area product of the MTJ considerably compared to ours. These two factors might be the 
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cause for the 650 times lower dissipation figure reported in Ref. [8] compared to what 
we find. Even then, the lower dissipation reported in Ref. [8] is still 45 times higher 
than that encountered in the straintronic spin neuron. If we carry out the comparison 
between similar designs (in-plane anisotropy magnets, similar energy barrier heights to 
maintain similar resilience to thermally induced random firing), then the difference is 
even more stark; the straintronic spin neuron is 29445 times more energy-efficient and 
yet 10 times faster. The current-driven spin neuron however has however one small 
advantage; it has a 50% higher on/off ratio of the MTJ resistance because the angular 
separation between the two stable orientations of the soft magnet in the MTJ is ~ 180° 
for the current-driven spin neuron and ~ 90° for the straintronic spin neuron. 
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-1-'-1-1-1-i-1->-1- 1 

0.00 0.05 0.10 0.15 0.20 

Voltage (V) 

Figure 5. Transfer characteristic (or firing behavior) for the straintronic spin neuron 
at room temperature in the presence of thermal noise. Results are shown for positive 
threshold only since no additional information can be gleaned from the negative 
threshold segment. Since the simulation is terminated immediately upon completion 
of firing, no fluctuations in the transfer characteristic are visible in the low-resistance 
state. 



5. Room-temperature firing behavior of spin neurons 

The room temperature firing characteristics are found by solving the stochastic LLG 
equation which will generate a distribution of characteristics since each one is slightly 
different in the presence of thermal noise. Each characteristic will show slightly different 
switching threshold and slightly different switching delay. 

Figure 5 shows the R(V)/Rl versus V characteristic of the straintronic spin neuron 
at 300 K when switching from high- to low-resistance state (the characteristic for 
switching from low- to high-resistance state will be qualitatively similar and hence not 
shown). As many as 10,000 switching trajectories were simulated from the stochastic 
LLG equation and make up this plot. There average switching time was 0.61 ns when the 
voltage appearing at node P in Figure 1 was above the 0 K threshold voltage of 93.15 mV. 
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Figure 6. Transfer characteristic (or firing behavior) for the current-driven spin 
neuron at room temperature in the presence of thermal noise. 


Clearly, the threshold has been significantly broadened at room temperature, indicating 
that there is a significant probability of pre-mature firing (firing before reaching the 
threshold defined at OK) because of random thermal torque, as well as failure to fire at 
or slightly beyond the threshold for the same reason. The likelihood of false firing is a 
matter of concern which calls for further investigation of thermal degradation. 

One measure of the threshold broadening is the ratio AV/Vth, where is the 
average voltage at which the neuron fires and AV is the standard deviation. For the 
straintronic spin neuron studied, this quantity turns out to be 16.5%. Of course this can 
be reduced by increasing Vth (by making the soft nanomagnet more shape-anisotropic 
to increase the in-plane shape anisotropy energy), but this will also increase the energy 
dissipation which varies roughly as V t 2 h (because the dissipation is dominated by the CV 2 
loss). Let us say that 1% broadening is acceptable. Then we will have to increase 
16.5-fold (assuming that it does not change AH), resulting in an energy dissipation of 
8.83aJ x 16.5 2 = 2.4 fj. This is still much less than what is encountered in CMOS-base 
implementations (0.7 pJ reported in Ref. [8]). 

The same thermal broadening is not only present in a current-driven spin neuron, 
but it is much worse there. Figure 6 shows R(I)/Rl versus / for the current-driven spin 
neuron where I is the driving current. Here, 5,000 switching trajectories were simulated 
to produce this plot because simulation of 10,000 trajectories became computationally 
prohibitive (because of the longer simulation time). The average switching time was 
now 7.26 ns above the 0 K threshold current of 127.13 /i A. The quantity A I / Ith 
turns out to be 69.92% for the simulated neuron. Had we been able to simulate more 
trajectories, we might have found the broadening to be even worse. Once again, the 
relative broadening can be reduced by increasing I t h (by increasing the shape anisotropy 
of the soft nanomagnet) but of course at the cost of increasing dissipation since the 
latter varies roughly as I 2 h (because the dissipation is dominated by the I 2 R loss). Once 
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again, if 1% broadening is acceptable, then we will need to increase the threshold 69-fold 
(assuming this does not change A/), resulting in an energy dissipation of 0.26pJ x 69 2 = 
1.23 nJ. This would make it 1768 times worse than CMOS-based neurons and therefore 
the viability of STT-based spin neurons at room temperature is questionable. Even if 
better design (use of spin Hall effect, perpendicular anisotropy magnets, etc.) can reduce 
the energy dissipation by two orders of magnitude, it will still provide little advantage 
over CMOS implementations. In contrast, the straintronic spin neuron is still ~290 
times more energy-efficient than its CMOS-based counterpart at room temperature. 

6. Conclusion 

In conclusion, we have proposed and analyzed a straintronic spin-neuron that is orders 
of magnitude more energy-efficient than the usual current-driven spin neuron and also 
faster. The primary obstacle they both face is the significant broadening of the firing 
threshold at room temperature. For the current-driven spin neuron, it may not be 
possible to mitigate this problem without increasing the energy dissipation to the point 
where it is no longer superior to CMOS-based implementations. Fortunately, for the 
straintronic spin neuron, it may still be possible to mitigate this problem while retaining 
a significant energy advantage over CMOS. 

7. Appendix 

The broadening in Fig. 6 for current-driven spin neuron due to thermal noise at room 
temperature can be reduced if we choose a different soft material with larger Gilbert 
damping coefficient (a). Since the switching current density is approximately linearly 
proportional to a, the energy dissipation will increase quadratically. In order to illustrate 
this, we chose a hypothetical material which is identical to CoFeB in all respects, except 
its Gilbert damping coefficient a is 0.1 and carried out the stochastic LLG simulations. 
The switching current turned out to be 1.44 mA for a current pulse duration of 10 
ns, but the broadening (AI/I^) was reduced to 11.25%. However, the extremely large 
switching current results in energy dissipation of 23.8 pJ which is clearly prohibitive. 
This case study shows us that although thermal degradation may be countered with 
increased energy dissipation, the price to be paid in energy may be prohibitive. 
Programmable synapses: In this paper, our synapses have fixed weights because 
they are implemented with passive resistors. Programmable synapses are required for 
more versatile renditions of neural computing such as Spike Time Dependent Plasticity 
(STDP) models of neural networks [50, 51, 52, 53] that are popular for Hebbian learning 
[54], For STDP models, the synapse weight should be programmed by a spike signal, 
e.g. it should be a function of a time-integrated spike current. The time integrated 
spike current is a charge that could be stored in a capacitor which then applies a 
voltage across a piezoelectric layer that generates strain and rotates the magnetization 
of a magnetostrictive magnet that is elastically coupled to the piezoelectric layer. The 
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magnetostrictive magnet could be the soft layer of a magneto-tunneling junction whose 
resistance is thus programmed by the time integrated spike signal, resulting in the 
appropriate programmable synapse. This is a different topic and would be treated 
elsewhere. 
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